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Abstract: We consider a real-time multi-ierver system with homogeneous servers (such as overhearing devices, 
unmanned aerial vehicles, machine controllers, etc.) which can be maintained/programmed for different kinds 
of activities (e.g. passive or active). This system provides a sendee for real-time tasks arriving via several 
channels (such as communication channels, surveillance regions, assembly lines, etc.) and involves 
maintenance. We address the worst case analysis of the system working under maximum load with preemptive 
priorities assigned for senders of different activity type. We consider a model with ample maintenance facilities 
and single joint queue to all channels. 

We provide various analytical approximations of steady state probabilities for these real-time systems, discuss 
their quality, compare the results and choose the best one. 
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I. Introduction 

Real-time systems (RTS) are imbedded in most modern technological structures, such as self-guided 
missiles, production control systems, aircraft and space stations, reconnaissance, radars, robotic and 
telecommunications systems, etc. 

According to [ 1 ]: “real time systems are those systems that produce results in a timely manner”, i.e. an 
action performed out of time limits (too late or too early) may be useless, and sometimes harmful - even if such 
an action or computation is functionally correct. 

■ a . J , -■ 

The data arriving to RTS is mostly incomplete or uncertain and the output may also be uncertain. That is 
why many of RTS include stochastic as well as dynamic components. 

Many scientific communities were treating various RTS models, and there exists a rich literature covering 
this area. 

We will focus on RTS with a zero deadline for the beginning of job processing. In these systems, jobs are 
executed immediately upon arrival, conditional on system availability. That part of the job which is not 
executed immediately is lost forever, and queueing of jobs (or their parts) in such systems is impossible. 

The following works treat this kind of RTS. In ([2], [3]) it was proven that the non-mix policy of never 
relieving an operative server maximizes the availability of a multiserver single-channel RTS involving 
preventive maintenance and working in general regime with any arrival pattern under consideration and constant 
service and maintenance times. In [4] and [5] multi-server and multi-channel (identical servers and channels) 
RTS (with unrestricted and restricted number of maintenance facilities respectively), working under maximum 
load regime were treated as finite source queues ([6]). In [7] various performance measures for RTS with 
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arbitrary number of channels operating under a maximum load regime were presented. In [8] and [9] multi¬ 
server and multi-channel RTS working in general regime were studied. 

In [10] it was shown that even very large number of servers in RTS with ample maintenance facilities does 
not guarantee the maximum system availability, and optimal routing probabilities were computed analytically 
(for exponentially distributed service times) and via Cross Entropy (CE) [11],[12], [13] simulation approach 
(for generally distributed service times). These results were extended for RTS with limited maintenance 
facilities in [14]. 

RTS with priorities were studied in [15], [16] (preemptive) and [17] (non-preemptive) respectively. 

In [16] a multi-server and multi-channel RTS with single joint queue of servers for each channel and 
preemptive priorities was studied and a set of balance linear equations for steady-state probabilities was 
obtained. Unfortunately, these equations do not have analytical solutions. 

The work presented here provides various analytical approximations of steady state probabilities, using the 
modifications of techniques proposed in [18] and [19]. We compare the results and choose the best one. 

/; \\ 

The paper is organized as follows: In Section 2, the description of the model with preemptive priorities is 
presented. Section 3 provides balance equations for model with ample maintenance teams. Section 4 several 
approximation methods are described. In Section 5 some numerical results are presented. Finally, Section 6 is 
devoted to conclusions. 


II. The model 

The most important characteristics of RTS with a zero deadline for the beginning of job processing are 
summarized in [8]. A real-world problem was studied in [16]. Here we provide the formal description of RTS 
from [17]. 

The system consists of r identical channels. For proper performance each channel needs exactly one fixed 
server at any instant (maximum load), otherwise the information in this channel (at this specific moment) is 
lost. There are N servers (which are subject to breakdowns) in the system. A server, which is out of order, needs 
time units of maintenance. After repair a fixed server may be of u-th type/class of quality with probability 

p (u=l, These probabilities can be used as control parameters. Only after the repair is completed, the 

quality control procedure determines the quality type of fixed server. The fixed server of u-th type is operative 
for a period of time S u before requiring R i hours of repair. S u and R i are independent exponentially 

distributed random variables with parameters JLl u (u=l, and A respectively. It is assumed that there arc K 


identical maintenance facilities/teams in the system. Each team can repair exactly one server at a time. We 
consider the case of ample maintenance facilities (K > N), so that all N servers can be repaired 
simultaneously, if necessary. The duration times R ; of repair are i.r.v. exponentially distributed with parameter 

A , which does not depend on the quality type of the server (neither before nor after the repair). After repair, the 
fixed server will either be on stand-by or operating inside the channel. There is a single joint queue of fixed 
servers to all r channels. 

We assume that servers of the first kind of quality type have the highest priority, servers of the second 
quality type are the next priority in line, and so on. Finally, servers of the m-th quality type have the lowest 
priority. Server operating in any channel is interrupted, if another fixed server of higher priority type arrives 
from maintenance. When the operating server must be repaired, the fixed server with highest priority takes its 
place. 

The system works under a maximum load (worst case) of nonstop data arrival, which is equivalent to the case 
of a unique job of infinite duration in each channel (a total of exactly r jobs in the whole system). Thus, the 
nonstop operation of the channel is needed. 

If, during some period of time of length T, there is no fixed server to provide the proper operation of the 
channel, we will say that the part of the job /information of length T is lost forever. 
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III. Balance equations for steady state probabilities 

In [16] the state of the system was defined as {ti l n m ) , where n u , u=l,...,m is a number of fixed 

m 

servers of n-th quality type (obviously, ^ n u <N and n u > 0 ), and the corresponding steady state 


U= 1 


probability was denoted n j. There are 


N + m' 


v m 


states in total. The rate of arrival of fixed servers of u- 


th quality type from maintenance is denoted X u — /.p u , . This system can be represented as a network with two 
nodes. We will focus on the case r < N , otherwise all fixed servers will be busy and preemptive priority 
regime, therefore, will not work. 

The model with ample maintenance teams (K > N) was studied, and the following set of linear equations 
for steady state probabilities was obtained in [16]. 
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IV. Approximations 


fC?, 


In this section we study several approximations ([18], [19]) for the model under consideration. These 
approximations instead of Global Balance Equations (1) use Local Balance Equations (as is shown further), and 
provide analytical product form solutions [6] for steady state probabilities. We discuss their quality and choose 
the best one. We will use an RSS (root of sum of squares) criteria in comparison between different 
approximations: 


RSS 


N N-n , £ 

I z x - x (**o 

| 77 1 =0 n 2 =0 n m =0 




( 2 ) 


where P^ fl ( and P(n t ,..., n m ) are exact numerical solution (of equations (1)) and analytical approximation 
respectively. 
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We will use the state transition rate diagram on Fig. 1 (see Fig 2, [16]) in order to explain approximation 
methods. 


Fig. 1. State transition diagram of the system with m = 2, r = 1 < N, K > N 
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Method 1 

According to [19], in order to arrive from the state (/7|, n^) to the state (n\,n2) , where (/7j <n\) and 

(JI 2 — n 2 ), we move first along the horizontal line /7] to the state 2 ) using the local vertical cuts 

between neighbouring states for local equations. 

Then we move along the vertical column n 2 to the state (n\,H2) using the local horizontal cuts between 
neighbouring states for local equations. 
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To be more specific, we will show the route between the states (0,0) and (2,2): 

(0,0)->-(0,1)->(0,2)->(1,2)->( 2, 2) as it is shown in Figure 2, which is a part of Fig. 1. 


Figure 2. The route between the states (0,0) and (2,2) 
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The corresponding local balance equations between neighbouring states will be as follows: 
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Finally, using the mathematical induction, we obtain 
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It can be easily seen that only two-way arcs are used in this process, while one-way arcs ( like those between 
(1,0) and (1,1)) do not participate at all. But namely one-way arcs express the preemptive priority regime in the 
system. Therefore this approximation method must be adjusted for the system under consideration. 


Method 2 , 

This method is a modification of Method 1. The idea is to arrive first to some of the neighbouring states of 
(n \ ,112) via two-way arcs, and only then to use one-way arcs for local balance equation for (n \, U2) ■ 

In our example to reach from the state (0,0) to the state (2,2) it will work as follows. First to arrive to the states 
(1,2) and (2,1) according to Method 1: 
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(0,0) ->(0,1) ->(0,2) ->(1,2) 
(0,0)->(0,l)-Kl,l)->(2,l) 
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Next considering the state (2,2) and its links with the states (1,2) and (2,1) only, we obtain 
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Method 3. 

This method is a combination of the methods presented in [18] and [19] .First we are treating the states on the 
upper line // ( =0, i.e. (0,0), (0,1), ..., (0,N), while using only vertical cuts for local balance equations. Then 
we step down to the next horizontal line Pj = 1, while using diagonal cuts for local balance equations (see Fig. 

3). And so on we continue tovards the last horizontal line = N . 

The corresponding local balance equations are as follows: 
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for the line Pj = 2. 

Thus we arrive to the state (2,2) via the following route 

( 0 , 0 )->( 0 , 1 )->( 0 , 2 )->( 1 , 2 )->( 2 , 2 ). 

Solving these equations sequentially, we get: 
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Figure 3. Local cuts in the system with m = 2, r = 1 < N, K > N 
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for Method 3 correspondingly. 
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V. Numerical results 

In this Section, we present some numerical results, which allow us to compare between approximation 
methods 1-3 and to find the best one in terms of RSS. 

Table 1 contains RSS in the case of ample maintenance facilities and r = 1; N = 4,5,6; 1 = 10; m = 2; jui = 8, ji 2 
= 6; p, = 0.6, p 2 = 0.4. 

Table 1. Comparison between three approximation methods w.r.t. RSS. 


N 

Method 1 

Method 2 

Method 3 
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0.297 

0.299 

0:157 

5 

0.336 

0.334 

0.163 

6 

0.336 

0.344 

0.190 


' \ i 

It can be easily seen from the Table 1, that Method 3 has a minimal value of RSS, and is therefore the best one. 

Tables 2-4 contain values of steady state probabilities of all three approximation methods as well as exact results 
in the case of ample maintenance facilities and r = 1; IV = 4,5,6; ). = 10; m = 2; /(| = 8, p 2 = 6; p\ = 0.6, p 2 = 0.4. 


Table 2. Steady state probabilities obtained by approximation methods and exact solution, N = 4 
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Table 3. Steady state probabilities obtained by approximation methods and exact solution, N = 5 


Exact 

solution 

Method 

3 

Method 

2 

Method 

1 

States 

(n,,n 2 ) 

Exact 

solution 

Method 

3 

Method 

2 

Method 

1 

States 

(n h n 2 ) 

0.003 

0.020 

0.027 

0.025 
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0.005 

0.002 

0.002 

0,0 

0.024 

0.050 
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0.049 

1,2 

0.004 

0.006 

0.002 

0.008 

0,1 

0.121 

0.088 
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0.066 

1,3 

0.023 

0.066 

0.021 

0.022 

0,2 

0.217 

0.088 

0.048 

0.044 
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0.096 

0.132 

0.041 

0.044 
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0.209 
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0.055 
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0.000 
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0.009 
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0.007 
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0.000 

0.008 
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0.069 

0.062 

4,1 

0.000 

0.008 

0.079 

0.094 

4,0 


0.000 

0.006 

0.059 

0.070 
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Table 4. Steady state probabilities obtained by approximation methods and exact solution, N = 6 
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Method 
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Method 
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Method 
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0.003 

0.022 

0.021 

0.018 
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0.000 

0.001 

0.000 

0.000 
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0.084 
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0.001 
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0.004 

0.018 

0.004 
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0.011 

0.016 

0,3 

0.000 

0.008 

0.024 

0.02 

2,1 

0.096 

0.097 

0.021 

0.032 

0,4 

0.003 

0.025 

0.047 

0.04 

2,2 

0.209 

0.129 

0.028 

0.042 

0,5 

0.022 

0.041 

0.047 

0.04 

2,3 

0.139 

0.086 

0.019 

0.028 

0,6 

0.104 

0.073 

0.042 

0.036 

2,4 

0.000 

0.002 

0.001 

0.002 

1,0 

0.000 

0.009 

0.053 

0.045 

3,1 

0.000 

0.002 

0.005 

0.008 

2,0 

0.004 

0.054 

0.063 

0.054 

3,2 

0.000 

0.003 

0.015 

0.023 

3,0 

0.022 

0.025 

0.071 

0.06 

3,3 

0.000 

0.003 

0.034 

0.051 

4,0 

0.000 

0.009 

0.08 

0.068 

4,1 

0.000 

0.003 

0.051 

0.077 

5,0 

0.003 

0.018 

0.053 

0.045 

4,2 

0.000 

0.002 

0.038 

0.057 

6,0 

0.000 

0.007 

0.06 

0.051 

5,1 

0.000 

0.007 

0.008 

0.007 

1,1 


VI. Conclusions 

We have found a good analytical approximations (Method 3, equations (5) and (8)), for RTS with a single 
joint queue, ample maintenance facilities and preemptive priorities working under maximal load. In our future 
research we will try to find good analytical approximations for RTS with a single joint queue and shortage of 
maintenance teams and for RTS with separate queue for each channel. 
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